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ABSTRACT 

I show that for a substantial fraction of planets detected in a space-based 
survey, it would be possible to measure the planet and host masses and distances, 
if the survey satellite were placed in geosynchronous orbit. Such an orbit would 
enable measurement of the microlens parallax tte for events with moderately 
low impact parameters, (3 < 0.05, which encompass a disproportionate share of 
planetary detections. Most planetary events yield a measurement of the angular 
Einstein radius ^e- Since the host mass is given by M = Oe/ktte where k is a 
constant, parallax measurements are the crucial missing link. I present simple 
analytic formulae that enable quick error estimates for observatories in circular 
orbits of arbitrary period and semi-major axis, and arbitrary orientation relative 
to the line of sight. The method requires photometric stability of ~ 10~^ on 
~ 1 orbit timescales. I show that the satellite data themselves can provide a 
rigorous test of the accuracy of the parallax measurements in the face of unknown 
systematics and stellar variability, even at this extreme level. 

Subject headings: gravitational lensing: micro — planetary systems 



1. Introduction 

Microlens planet detections routinely yield the planet-star mass q = m/M, but since the 
host star is generally very difficult to observe, the host mass M (and so the planet mass m) 
and the system distance Di usually remain undetermined. In principle, one can determine 
these quantities from the microlensing event itself, provided th e angular Einstein radius 6-^ 



and the microlens parallax tte are both measured (lGouldlll992l ) 



, . 6'e „ 46* o 1 Kias , 

where tte = AU/fE is the ratio of the Earth's orbit to the Einstein radius projected on 
the observer plane, and tTi-ci = AU[Dj^^ — D^^] is the lens-source relative parallax, with Ds 
generally being known from direct observation of the source. 
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To date, 6^ has been measured in the majority of planetary events, even though it is 
rarely (~ 0.1%) measured in ordinary events. This is because planets are detected when the 
source passes near or over a caustic induced by the planet, which yields a measurement of 
the source size relative to the Einstein radius. 

Hence, if a method could be found to measure tte for a large fraction of events, the 
information returned from microlensing planet searches would be radically improved. The 
microlens parallax is actually a vector tte, with the direction being that of the lens-so urce 
relative motion. It can in principle be measured either from the gro und (jGouldlll992l ). by 
combining spa ce-based and ground-based observations (lRefsdallll966l ). or from space alone 
dHonmalflOOol ). 



When I proposed the pure ground-based method (lGouldlll992l ). I illustrated it with an 
event that lasted 4 years, so that the lightcurve oscillated annually around the standard mi- 
crolensing form. In practice, such 4- year events are never seen, but tte can be measured from 
less regular features in shorter events. However, typical events have timesca l es ~ 20 days, 
so that the great majority are too short to permit a measurement. iHonmal (119991 ) proposed 
a purely space-based method using Hubble Space Telescope (HST) observations. Since the 
HST orbital period is only 1.5 hours, the oscillations are obviously much shorter than ^e- 
However, the HST orbital radius is only i?®, whereas typically te ~ 1 — 10 AU. Hence, 



the oscillations would typically be minuscule. To overcome this obstacle, iHonmal (119991 ) 
proposed that the observations should be made during caustic crossings, which he showed 
dramatically increase the strength of the signal. However, alerting HST to imminent caustic 
crossings is extremely difficult, more so for planetary events, and in fact such measurements 
have never b een carr i ed ou t, nor even attempted. The ground-space combination originally 
proposed by iRefsdall (119661) remains fe asible in principle, but each event must be observed 
very frequent l y (jGaudi fc Gouldl 119971 ) . implying a huge space-craft investment. Recently, 
Gould fc Yed (120121 ) showed that a much simpler and cheaper ground-space approach was 
feasible for high-magnification events, but these constitute a small subset. 

Here I show that a single survey satellite in geosynchronous orbit is a near-optimal 
solution to this problem. As in my original proposal, oscillations are induced on a stan- 
dard microlensing lightcurve, but since these have periods of one day rather than one 
year, they are sho rter tha n the e vent. The amplitude of the oscillations is 6.6 times larger 
than those of the iHonmal (119991 ) proposed HST observations. Moreover, no special warn- 
ing system is required. It is true that this approach only works for moderately high- 
magnifi cation {A > 10) eve nts, but these are more sensitive to planets than more typical 
events (iGould &: Loeblll992l ). Hence, this approach can yield parallaxes for a large fraction 
of planetary events. Moreover, such a geosynchronous survey satellite is a very real possi- 
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bility. The iDecadal Survey Committed (|2010[ ) report recommended a Wide Field Infrared 
Space Telescope (WFIRST) as its highest priority. A large fraction of the observing time 
would be devoted to microlensing planet searches. Although originally recommended for L2 
orbit, a geosynchronous orbit is now being actively considered. Hence, it is most timely to 
investigate the implications of such an orbit for microlensing parallaxes. 



2. Analytic Treatment of Single-Observatory Parallax 

Consider an observatory in circular orbit with period P and semi-major axis a. If this is 
an Earth orbit, these two quantities are obviously related, but here I begin with the general 
case. Let the microlensing target be at a latitude A with respect to the plane of this orbit, 
so that the projected semi-minor axis is 6 = a sin A. I then normalize these axes to an AU 

a a sin A 

^""^=AU' ^^ = ^' 

Now consider that continuous observations of a microlensing event are made at a rate 
per Einstein timescale ts, with signal-limited flux errors 

a' = alA. (3) 

This condition is not as restrictive as may first appear. We are only concerned with scaling 
of the errors when the source is relatively highly magnified, so that this form need not hold 
all the way down to baseline, where sky and blending may be important. Let the lens pass 
the source on its right at an angle Q relative to the major projected axis of the satellite orbit, 
with impact parameter /3 (in units of ^e)- Then the lens-source separation vector u (in units 
of 6'e) is given by 

u = (r cos Q — (i ^m.Q ^ eyTrE cos(ci;t + 0), r sin 6' + /3 cos 6' + e^vTE sin(a;t + 0)) (4) 

where r = it — t^jt-^, to is the time of closest approach, uj = 27r/P, and (j) is the orbit phase 
relative to peak. By definition, tte = (vteji, tte.j.) = 7rE(cos ^, sin 6*). Hence, 

du e\\T cos{ojt + (j)) + e±/3 sm{ojt + (j)) du —e\\f3cos{cot + (f)) + e±Tsm{ut + (j)) 



Otie II u diiE ± u 

(5) 



The magnification is given by (lEinsteinl 119361 ) 



uV^FTl' du m(«2 + 2)(u2 + 4)' ^ ' 
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so that the derivatives of the flux with respect to the parallax parameters are 

dF ^dlnA ^ du dF ^d\nA ^ du 

= ' ^ = ' ^ 

where F^ is the unmagnified source flux. In order to make a Fisher-matrix analysis, I first 
assume that the magnification does not change much over an orbit. This will ultimately 
restrict the result to events for which tcs = PtE > P- I will discuss the more general case 
further below. I then evaluate the mean contribution to the Fisher matrix, averaged over 
one orbit 

dF \^ F^ ,/ du \\ ^/d\nA\^ F^ 6A{ey + elP^) 



1 OF dF F^ 64(ei-e^)r/3 

a^^dTTE,l\dTTE,±' m5(m2+ 4)5/2(^2 + 2) ^ ' 

and so obtain the inverse covariance matrix by integrating over all observations 



NFl f elG2 + elGo 



^ = ^( " n , ) (11) 

where 



roo 64 X"/?^"" 

Gn{f3)=J^ dx ^ ^2)5/2(^2 + ^2 + 4)5/2(^2 + /32 + 2) ' ^^^^ 

In the limit /3 <^ 1, it is straightforward to show that 

Go(/3)^^r', G2(/3)-^ir'. (13) 

Figure [1] shows that this approximation holds quite well for /3 < 0.1 (i.e., A^ax ^ 10). Using 
this approximation, we obtain. 



/Too^ /(l + 2sin2A)-V2 
^(vrE,±); ViVF, e V (2 + sin^ A)-!/^ 



(14) 



Since the ratio of the quantities in the right-hand vector is never greater than \/2 (and is 
■\/3/2 for A = 30°, which is typical of Galactic bulge fields and an equatorial orbit), the 
parallax error is reasonably well characterized by the harmonic rms of these two values, 

. (£ii±ai)-"^ = o.o23( J^L)-"Vi±£E!A)-"V V^k) ( JL) 

^ ^ V 2 / VlOOOO/ V 1.25 / V 0.01 / V0.05/ V6.6i?e/AU 

(15) 
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3. Application to WFIRST 

The fiducial values in Equation f lTSj) are plausible for a WFIRST style mission. For 
example, a continuous sequence of 3-minute exposures yields N ~ 10000 for a = 20 day 
event. A large fraction of sources would have 1% errors in such an exposure, whether 
continuous or co-added. Note that the formula does not really specify ctq and N separately, 
but just ao/y/N. Moreover, for typical = 20 day events and a P = 1 day orbit, the 
assumptions of the derivation basically apply for f3 = 0.05, i.e., tes = /3^e = P- Of course, the 
integral in Equation ( !T2|) goes to infinity, while the observations do not, but the overwhelming 
contribution to this integral is from a few tcs near peak. The fiducial parallax error in 
Equation ( fT5l) is quite adequate to measure the mass and distance of disk lenses but somewhat 
marginal for bulge lenses. For example, for M = 0.5 Mq, a typical disk lens at Di = 4kpc 
has tte ~ 0.18, but a typical bulge lens at Ds — = 0.75 kpc has tte ~ 0.06. At the very 
least, however, such a measurement would distinguish between bulge and disk lenses. 



4. Systematics and Stellar Variability 



Naively, Equation ( IT5l) appears to require "effective errors" ofN-^/^ao/Fs ~ 10-^ which 
must be achieved in the face of both astrophysical and instrumental effects. In particular, it 
would seem to challenge the systematics limit for real crowded-field photometry, even from 
space. However, what is actually required is not control of systematics to this level, but only 
control of systematic effects that correlate with orbital phase. This requirement is still not 
trivial for geosynchronous orbit because of phase-dependent temperature variations, but it 
is not as intractable as vetting all systematic errors at this level. 

Nothing is known about stellar variability at this level in H (likely WFIRST passband). 
Even Kepler has only a handful o f stars for which it can make measurements of this precision 
in its optical passband (Fig. 4 of iMcQuillan et al.l 120121 ). Ho wever, again what is rele vant is 
not variability per se, but variability on 1-day timescales, and iMcQuillan et al.l (|2012[ ) report 
that FGK stars typically vary on > 5-day timescales. Moreover, most forms of intri nsic 
variability are very subdued in H band relative to the optical. See iGould et aL J2013h for 
a spectacular example of this in a microlensing event. Finally, IMcQuillan et al.l ( 120121 ) find 



that variability declines with increasing proper motion (so presu mably age), and the bulge 
is mostly much older than the disk, though perhaps not entirely (IBensby et al.ll2013l ) 
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5. Rigorous Test of Parallax Accuracy 

One would nevertheless like a rigorous test that undetectable variability and systematic 
effects are no t corrupting the re sults. Fortunately the satellite data themselves provide 



such a test. iGould et al.l (Il994f ) showed that even very short microlensing events yield 
"one-dimensional parallaxes" from the asymmetry in the lightcurve induced by the Earth's 
(approximately constant) acceleration toward the Sun. That is, vrE,0,|| is well constrained 
while vrE,0,± is very poorly constrained. Normally, such 1-D parallaxes are not considered 
useful because the amplitude of vte (and so M) is also poorly constrained. However, in the 
present case these 1-D parallaxes serve two very important functions. Figure [2] shows the 
lightcurve of a simulated event with 72 days of observations centered on the Vernal Equinox 
in each of 5 years, similar to the anticipated schedule of WFIRST. The event is assumed to 
have the fiducial parameters from Equation f|T5l) : = 20 days, /3 = 0.05, N~^/'^{ao/ Fg) = 
10~^, with to right at the Vernal Equinox. The lower panel shows the difference between 
this event as observed from a geosynchronous orbit, and the same event observed from an 
inertial platform. Note that since the orbit is equatorial, (vteji, vrE,±) = (vrE,_B, 7rE,iv), i.e., 
the components of tte in the East and North directions projected on the sky. The event is 
at (a, 5) =(18:00:00,-30:00:00), and is assumed to have a parallax {ite,e,t^e,n) = (0.1,0.1). 
The oscillations near peak are due to the satellite's orbital motion, while the asymmetries 
in the wings are from the Earth's orbit. The latter are much larger, but yield only 1-D 
information. This is illustrated in Figure El which shows the error ellipses due to 1) Earth- 
only, 2) satellite-only, 3) Earth+satellite. Also shown is the result of the analytic calculation 
(which assumed observations extending to infinity) . Note that there are two sets of solutio ns 



in each case, corresponding to the (/3 — )■ —(3) degeneracy ( ISmith et al.ll2003l : lGouldll2004j ) 



The first point is that if 71e,e as derived from geosynchronous parallax agrees with the 
much more precise value derived from the Earth's orbit, then one can have good confidence 
in 7rE,Ar, for which there is no direct test. This is true on an event by event basis, but more 
true for the ensemble of parallax measurements. 

The second point is that these two parallax measurements will be automatically com- 
bined in the fit to any lightcurve, which means that tce^e will be much better determined than 
tte.at. Depending on the relative values of these two components, this may add important 
information in some cases. 

Finally, the ensemble of geosynchronous measurements of 7rE,E provides a test of the 
accuracy of the 1-D Earth-orbit measurement of this quantity from the wings of the event. 
This is important because the source stars may show variability on 10-day timescales, even 
if they do not vary on 1-day timescales. 
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Using techniques similar to those used to derive Equation f lMj) . it is straightforward 
to show that the 1-D parallax error due to Earth acceleration for short, relatively high- 
magnification events (and infinite observations) is 



where 58 days is one radian of the Earth's orbit and rj is the projected Earth-Sun separation 
in AU at time to- Since 77 ~ 1 for bulge observations made during the equinoxes, the ratio 
of geosynchronous-to-Earth parallax errors is 

0'(7rE)geosynch qQ V ^ - (17) 

cr(7rE,0,||)Earth V20day/ 0.05 1 day 1 day 

Since the effective pt-^ > P = 1 day, this implies that the Earth-orbit parallax will essentially 
always yield a precise check on the geosynchronous parallax in one direction. 



Discussion 



The derivation underlying Equation f lTSj) breaks down for t^s ^ P- the errors continue to 
decline with falling /3, but no longer linearly. They also become dependent on the orientation 
and phase of the orbit in a much more complicated way. From the present perspective, the 
main point is that the formula with /3 — P/^e provides an upper limit on the errors for 
events with yet higher peak magnification. 

Next, the errors derived here assume a point-lens event. However, since the observations 
would be near-continuous, it is likely that caustic crossings or near approaches would be 
captured. As pointed out by iHonmal (119991 ) such caustic effects can significantly enhance 
the signal. 

Another feature of these (and most) parallax measurements is that they work better at 
low mass, simply because tte oc M~^/^ is bigger. Space-based microlensing measurements 
have the potential to directly detect the lens when it is more massive (so, typically, brighter). 
For example, as the lens and source separate after (or before) the event, their joint light 
becomes extended and the centroids of the blue and red light sep arate (if the sourc e and l ens 
are different colors). These effects allowed iBennett et al.l (120061 ) and iDong et al.l (120091 ) to 
measure the host masses in two different planetary events using followup HST data. Because 
geosynchronous parallax works better at low mass, while photometric/astrometric methods 
work better at high mass, they are complementary. 



Finally, I note that such parallaxes would be of great interest in non-planetary events 
as well. Without 6e, such measurements do not yield masses and distances, but they do 
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serve as important inputs into Bayesian estimates of these quantities. Moreover, since the 
direction of tte is the same as that of the lens-source relative proper motion a parallax 
measurement provides an important constraint when trying to detect/measure the source- 
lens displacement away from the event. If the magnitude fi^ei can be measured from these 
data, then so can 9e = /irei^E, which in turn yields the mass and distance. 

This work was supported by NSF grant AST 1103471 and NASA grant NNX12AB99G. 
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Fig. 1. — Functional forms of Go{/3) and G2{/3) relative to the limiting forms given by 
Equation (IT^ . These approximations are excellent for (3 < 0.05 and very good for /3 < 0.1. 
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Fig. 2. — Simulated lightcurve of event with = 20 days, /3 — 0.05, to at the Vernal Equinox, 
and parallax 7rE,Ar = t^e,e = 0.1, observed from a geosynchronous equatorial orbit. The error 
bars are binned by day for display but the observations are assumed many times per day. 
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Fig. 3. — Error ellipses (1 a) for the event shown in Figure |2] using Earth-orbit-only (green), 
geosynchronous-only (red) and combined (black) information. The lightcurve asymmetry 
(Fig. [2]) due to the Earth's orbit yields only 1-D parallax information, but this serves as a 
critical check on the accuracy of the geosynchronous-orbit parallax. The secondary minimum 
at 7rE,iv ~ —0.07 is due to (/3 — )■ — /3) degeneracy. In this example, it is disfavored by 
Ax^ = 11 based on geosynchronous data and by 15 based on all data. The analytic error 
estimate (based on infinite data) is shown in cyan. 



